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(57) ABSTRACT 

A method for estimating the time varying mechanical prop- 
erties of a material comprising: propagating acoustic waves 
through a sample of material; measuring signals correspond- 
ing to acoustic waves after they propagate in the sample; and 
characterised by: comparing attributes from the measured 
signals with attributes from corresponding signals provided 
by a calculating model, the calculating model accounting for 
at least initially estimated acoustic properties of the material; 
using the differences between the measured and calculated 
signals to update the initially estimated acoustic properties 
of the material in the calculating model; processing N 
iterations of the comparison and updating either until the 
differences between measured and calculated signals are 
within a given tolerance factor or when the number i reaches 
a prescribed maximum; and calculating the time varying 
mechanical properties of the material from the final iterated 
acoustic properties. Also described is an apparatus for 
implementing the method. 
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ACOUSTIC METHOD FOR ESTIMATING 
MECHANICAL PROPERTIES OF A MATERIAL 
AND APPARATUS THEREFOR 

TECHNICAL FIELD 

[0001] The present invention relates to an acoustic method 
for estimating the mechanical properties of a material and an 
apparatus therefor. In particular, the invention provides a 
new method and apparatus for estimating the acoustic prop- 
erties of set cements used in an oil well or the like. 

BACKGROUND INFORMATION 

[0002] After drilling an oil well or the like, the annular 
space surrounding the casing is generally cemented. Such an 
operation is intended to consolidate the well, to protect the 
casing and also, essentially, to isolate geological layers so as 
to prevent fluid exchange between the various formation 
layers, where such exchange is made possible by the path 
formed by the drilled hole. The cementing operation is also 
intended to prevent gas from rising via the annular space and 
to limit the ingress of water into the production well. Good 
isolation is thus the primary objective of the majority of 
cementing operations carried out in oil wells or the like. 

[0003] Consequently, the selection of a cement formula- 
tion is a critical issue in such cementing operations. The 
appropriate cement formulation helps to achieve a durable 
zonal isolation, which in turn ensures a stable and productive 
well without requiring costly repair. Important parameters in 
assessing whether a cement formulation will be optimal for 
a particular well environment are the mechanical properties 
of the cement after it sets inside the annular region between 
casing and formation. Compressive and shear strengths 
constitute two important cement mechanical properties that 
can be related to the mechanical integrity of a cement sheath. 
These mechanical properties are related to the linear elastic 
parameters namely: Young's modulus, shear modulus, and 
Poisson's ratio. It is well known that these properties can be 
ascertained from knowledge of the cement density and the 
velocities of propagation of the compressional and shear 
acoustic waves inside said cement. Thus, ultrasonic mea- 
surements can be used: these non invasive measurements are 
particularly interesting as opposed to other well known 
mechanical techniques whereby samples are stressed to a 
failure stage to determine their compressive or shear 
strength. 

[0004] Acoustic tools are usually used to perform these 
acoustic measurements. These tools are lowered inside a 
well to evaluate the cement integrity through the casing. 
However, interpretation of the acquired data is particularly 
difficult, especially due to the various compositions of the 
cements and data derivation induced by — among other — 
electronics inside said acoustic tools. Consequently, several 
mathematical models have already been developed to simu- 
late the measurements, those models being very helpful to 
anticipate the performance of the evaluation tools as well as 
to interpret the tool data. However, those models necessitate 
the knowledge of both the velocities and attenuations of 
acoustic waves that propagate inside the cement. Thus, 
providing accurate values to those models is a very impor- 
tant aspect in interpreting the measurements of the different 
acoustic tools. 

[0005] U.S. Pat. No. 4,259,868, to Rao et al, describes a 
device that measures compressive ultrasonic waves transit 



time through cement samples inserted in an auto-clave held 
at constant temperature and pressure. The transit time infor- 
mation is converted to compressive strength on the basis of 
established relationships between the two quantities. These 
relationships are obtained empirically from carrying out 
ultrasonic transit time measurements and mechanical com- 
pressive strength measurements on many samples. As far as 
the transit time of compressive ultrasonic waves decreases 
along the setting of the cement, the device, performed on a 
cement slurry, provides the cement set time, or the time it 
takes for the slurry to gain a minimal compressive strength, 
and the time for development of full compressive strength. 
These parameters are used as inputs for efficient scheduling 
and carrying out of a well cementing job. However, the 
reliability and accuracy of this device to provide cement 
compressive strength relies on the accuracy of the correla- 
tion relationships that relate compressive strength to com- 
pressional wave transit time. 

[0006] Another U.S. Pat. No. 5,741,971 to Lacy also 
provides the characteristics of cement slurries measured 
acoustically in an auto-clave. These characteristics include 
dynamic Young's modulus, density, static viscosity, com- 
pressive strength and expansion or contraction of set 
cement. The Young's modulus measurement is based on an 
estimation of a compressional wave transit time through a 
cement sample. The estimation of Young's modulus from 
the transit time information is based on what is believed to 
be an empirical relation between the two quantities. Fur- 
thermore, it is speculated that the derivation of this empirical 
relation is based on calibration measurements for cements 
with well known or independently measured Young's 
moduli. This patent also describes an additional empirical 
relation to relate the estimated Young's modulus to the 
unconfined compressive strength of the tested cement. This 
relation is derived using best fit approaches applied to data 
relating the two quantities where the compressive strength is 
measured independently using a conventional measurement. 

[0007] All of these embodiments of ultrasonic measure- 
ments are geared towards propagating and detecting com- 
pressional waves only. However, since shear waves only 
propagate through solid media, detecting shear waves in the 
slurry as a function of time indicates when the transition 
from liquid to solid/gel occurs, which is an important input 
in determining the mechanical properties of a cement. In 
U.S. Pat. No. 5,412,990, granted to D'Angelo et al, the onset 
of shear waves transmitted through a cement slurry is 
monitored and directly related to the cement thickening 
time. Nevertheless, as this method does not take compres- 
sional waves into account, it suffers from being inaccurate. 

SUMMARY OF THE INVENTION 

[0008] A principle aim of the invention is to propose a 
method suitable for estimating the time varying mechanical 
properties of a material, said method being particularly 
accurate at providing reliable values for said time varying 
mechanical properties. 

[0009] To this is end, the invention provides a method for 
estimating the time varying mechanical properties of a 
material comprising: 

[0010] propagating acoustic waves through a sample 
of said material; 
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[0011] measuring signals corresponding to said 
acoustic waves after they have propagated in said 
sample. According to the invention, the method 
comprises: 

[0012] comparing attributes from said measured 
signals with attributes from corresponding signals 
provided by a calculating model, said calculating 
model accounting for at least initially estimated 
acoustic properties of said material; 

[0013] using the differences between said mea- 
sured and calculated signals to update said initially 
estimated acoustic properties of said material in 
the calculating model; 

[0014] processing (N) iterations of said compari- 
sons and updatings until a given calculating stop; 

[0015] calculating the time varying mechanical 
properties of said material from said final iterated 
acoustic properties. 

[0016] The method according to the invention provides 
particularly accurate values for the time varying mechanical 
properties of the material because it is based on a iteration 
process towards the measured values, said iteration process 
being based upon simple and reliable parameters. 

[0017] In a preferred embodiment, the acoustic waves are 
both compressional and shear waves, which emphasises the 
efficiency of the method because the correlation of the 
different values permits the calculated mechanical properties 
of the material to be checked. 

[0018] According to an example of the method, the mea- 
sured and calculated signals are voltage values as a function 
of time and the attributes of the measured and calculated 
signals are the transit times of the acoustic waves and 
amplitudes of said signals. Consequently, the measured 
signals are very easy to acquire with conventional acoustic 
tools. 

[0019] The acoustic properties of the material in a pre- 
ferred embodiment of the method comprise the velocity of 
the acoustic waves through the sample of said material. 
These properties thus allow the easy calculation of different 
mechanical properties of the material such as the Young's 
modulus, Poisson's ratio and shear modulus. 

[0020] In a preferred embodiment of the method, at the 1 th 
iteration (i=N), the value of the velocity for the acoustic 
waves through the sample of material is updated according 
to: 



V/ +l = Vj ■ where S Y . = --^ 5 V* 



[0021] where 1 means either the shear (S) or compressional 
(P) waves, t" f is the transit time of the waves in the 
measured signals, tj is the corresponding transit time of the 
calculated signal with a certain value V, for the wave 
velocity, d is the distance between the propagating and the 
measuring means, h is the sample thickness, and V f is the 
sound speed in the medium. 



[0022] According to this preferred embodiment, the itera- 
tion stops when the relative difference between two succes- 
sive iterated values of the velocity for the acoustic waves 
through the sample of material is within a prescribed toler- 
ance factor e v , such as: 




[0023] In another preferred embodiment of the method of 
the invention, the acoustic properties of the material com- 
prise attenuation of the acoustic waves through the sample 
of material. This permits interpretation of data acquired by 
acoustic tools in in situ conditions, such as down an oil well. 

[0024] In this preferred embodiment, at the 1 th iteration 
(i=N)» the value of the attenuation for the acoustic wave is 
updated according to: 



or*/" 1 = aj + 6 Ql ; where 6 Ql = — 1 ^ ' a*, 



[0025] where 1 means either the shear (S) or compressional 
(P) waves, A/** is the amplitude of the acoustic waves in the 
measured signals and A x is the corresponding amplitude of 
the calculated signal with a certain value for the attenuation 
a,. 

[0026] According to this embodiment, the iteration stops 
when the relative difference between two successive iterated 
values of the attenuation for the acoustic wave is within a 
prescribed tolerance factor e ttJ such as: 



[0027] In a preferred embodiment of the method of the 
invention, the propagation and measuring means comprise 
ultrasonic compressional transducers, including a transmit- 
ter and a receiver. 

[0028] The use of non-contact transducers permits avoid- 
ance of shear coupling effects typically present when using 
contact transducers. Actually, the shear coupling between 
the transducers and the sample introduces in the measured 
signal amplitudes non -controllable variations that can't be 
accounted for in the calculating method. 

[0029] In another preferred embodiment of the method of 
the invention, the calculating model accounts for character- 
istics of the propagating and measuring means. These char- 
acteristics of the propagating and measuring means may 
comprise the radiation and the reception properties of said 
propagating and measuring means and the estimation of the 
transfer function of the electronic devices included in said 
means. The method according to the invention can also 
account for the effects of environmental parameters such as 
temperature and pressure. 

[0030] When accounting for the different characteristics of 
the propagating and measuring means, the method according 
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to the invention provides very reliable values of the 
mechanical characteristics of the tested material, even when 
environmental conditions are aggressive, for example under 
high temperature and/or high pressure, such as down an oil 
well. 

[0031] It is also an aim of the invention to provide an 
apparatus for propagating and measuring acoustic waves 
through a sample of material, the mechanical properties of 
which vary with time, said apparatus comprising: 

[0032] at least one transmitter (1) propagating said 
acoustic waves towards at least one receiver (2), said 
receiver detecting said acoustic waves after they 
have passed through said sample; and means for 
calculating properties of the sample based on the 
detected acoustic waves; characterised in that the 
apparatus comprises: 

[0033] means for varying the incidence angle of 
the acoustic waves on the sample of the material; 
and 

[0034] means for 

[0035] comparing attributes from said detected 
signals with attributes from corresponding sig- 
nals provided by a calculating model, said cal- 
culating model accounting for at least initially 
estimated acoustic properties of said material; 

[0036] using the differences between said mea- 
sured and calculated signals to update said 
initially estimated acoustic properties of said 
material in the calculating model; 

[0037] processing (N) iterations of said com- 
parisons and updatings until a given calculating 
stop; 

[0038] calculating the time varying mechanical 
properties of said material from said final iter- 
ated acoustic properties, 

[0039] This apparatus can provide the measured signals 
that will be compared to corresponding calculated signals 
according to the method of the invention. The fact that the 
incidence angle of the acoustic waves on the sample can be 
varied allows both compressional and shear waves to be 
propagated in the sample. This permits to provide very 
accurate values of the mechanical characteristics of the 
material to be provided when combined with the method 
according to the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0040] Further characteristics, advantages and details of 
the method and apparatus of the invention will be apparent 
from the description below which is made with reference to 
the accompanying drawings, given by way of example in 
which: 

[0041] FIG. 1 is a first embodiment of the measurement 
apparatus of the invention; 

[0042] FIG. 2 is a second embodiment of the measure- 
ment apparatus of the invention; 

[0043] FIG. 3 is a third embodiment of the measurement 
apparatus of the invention; 



[0044] FIG. 4 is a fourth embodiment of the measurement 
apparatus of the invention 

[0045] FIGS. 5fl and 5b are the structure and geometry of 
the embodiment represented in FIG. 1; 

[0046] FIG. 6 shows signals acquired with the embodi- 
ment of FIG. 1; 

[0047] FIGS. 7 to 9 show measured and model-calculated 
signals for compressional waves velocities and attenuation 
estimation; and 

[0048] FIGS. 10 to 12 show measured and model-calcu- 
lated data for shear waves velocities and attenuation signals. 

DETAILED DESCRIPTION 

[0049] As previously said, the present invention describes 
an acoustic method for estimating the time varying mechani- 
cal properties of a material. This method provides acoustic 
properties of a sample of said material, in order either to 
determine said time varying mechanical properties or to 
provide inputs in models anticipating the responses of 
acoustic evaluation tools, lowered for example in an oil well, 

[0050] The method consists in propagating acoustic waves 
through a sample of material to be tested. As represented on 
FIGS. 1 to 5, an apparatus allows said acoustic waves, 
namely both compressional and shear waves, to be propa- 
gated through a sample of material to be tested. The signals 
measured are then compared to corresponding signals pro- 
vided by a calculating model, said calculating model 
accounting for acoustic properties of the material and char- 
acteristics of the propagating and measuring means. The 
method then uses the differences between said measured and 
calculated signals to update the acoustic properties of said 
material in the calculating model and processes N iterations 
of said comparison and updating either until said differences 
between measured and calculated signals are within a given 
tolerance factor or when said number N reaches a prescribed 
maximum. Then, the method calculates the time varying 
mechanical properties of said material from said final iter- 
ated acoustic properties. 

[0051] In the first embodiment, illustrated in FIG. 1, a 
single ultrasonic transmitter 1 and a single ultrasonic 
receiver 2, both known in the state of the art (for example 
ultrasonically-pulsed acoustic piston transducers of about Va 
inch diameter), are facing and positioned at a distance d 
from each other. These transducers are surrounded by a 
medium 4 allowing the propagation of acoustic waves, said 
medium being water for example. A sample of material 3, 
some cement for example (or any material the mechanical 
properties of which vary with time), is placed between these 
two transducers. Preferably, these transducers are compres- 
sional ones (also known as immersion transducers) instead 
of contact transducers. Ultrasonic non-contact transducers 
are preferred because they do not suffer from shear coupling 
effects, occurring when using contact transducers. Actually, 
the shear coupling between the contact transducers and the 
sample introduces non-controllable variations in the signal 
amplitude and these variations cao significantly disturb 
mathematical models used to estimate such measurements. 
The coupling for non-contact transducers is only dependent 
on the material parameters of the fluid medium 4 and the 
cement sample 3, those parameters being accounted for in 
the latter described model. Using non-contact transducers 
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thus provide more accurate prediction of the propagation 
and detection of the acoustic waves through the sample. 

[0052] The distance d between the transmitter and the 
receiver is chosen so as to accommodate a sample 3 of 
material (for example between 10 mm and 20 mm). The 
thickness h of this sample is such that it can be placed 
between the transducers at stand-off s with respect to each 
transducer, in order to be able to rotate the sample or the 
transducers as it will be described below. The sample of 
material 3 can be a set cement sample having a rectangular 
shape with known or measurable thickness along which 
acoustic signals can propagate. In this case, the method 
according to the invention will provide data on the evolution 
of the mechanical characteristics of this sample when, for 
example, submitted to temperature or pressure. In an other 
embodiment, cement slurry may also be accommodated 
inside containers of rectangular shape. In this case, the 
method according to the invention allows to follow the 
evolution of mechanical characteristics of a cement that is 
thickening, and consequently the thickening time and the 
mechanical properties of this cement after setting. 

[0053] The cement sample is positioned preferably on a 
rotation stage such that its surfaces make an incidence angle 
theta (G) with respect to the transducers sensitivity lines as 
shown in FIG. 1. For given incidence angles, 8, of the 
cement sample, ultrasonic data is acquired. FIG. 6 shows 
typical signals acquired from a sample of class H cement 
(well known in the API system), for h=15 mm, d= 160.38 
mm and a slurry density p=2000 kg/m 3 . For small values of 
8 including zero, the signals consist of contributions from 
compressional wave propagation in the sample. As 8 
increases, contributions from shear wave propagation arise. 
The compressional and shear events are identified on the 
signals in FIG. 6. The amplitude and transit time of these 
events will then be used in the model-based method 
described below to extract the compressional and shear 
velocities and attenuations. 

[0054] FIG. 5a represents the structure of the apparatus 
according to the embodiment described in FIG. 1: a pair of 
piston transducers (preferably piezo-electric sensors) 
immersed in water is used in transmission mode to excite a 
layer of cement; data is acquired for normal incidence and 
oblique incidence 6 as the cement layer is rotated with 
respect to its mid -plane about a rotation axis located at equal 
distance between the transmitter and the receiver. FIG. 5b 
represents the geometry of the configuration modelled theo- 
retically: the transducers are aligned at an angle in the (x, z) 
plane. Their locations are 6-dependent. 

[0055] A second preferred embodiment of the apparatus, 
shown in FIG. 2, consists of using multiple transmitters 1 
and receivers 2 with a non-rotating sample 3. The transmit- 
ting and receiving elements are positioned and aligned at 
particular angles such that each transmitter-receiver pair has 
a beam incidence on the sample at a different angle. The 
number of pairs is chosen such that several incidence angles 
from zero to a maximum angle, 8 max, are covered, i.e. most 
of the acoustic waves are detected by the receivers after they 
have passed through the sample. Depending on the sample 
tested, 8 max may be near 40° for conventional cements and 
near 70° for cement with a large proportion of water. The 
large rotation angles are needed for this latter cements 
because their shear wave velocities are generally lower than 
the sound speed in the fluid. 



[0056] A third preferred embodiment of the apparatus, 
shown in FIG. 3, consists of a single transmitter 1 and a 
single receiver 2 with the cement sample 3 in a fixed 
position. The transmitter and receiver alignment angles 8 
may be controlled via an external mechanism. The signals at 
8=0 are used to determine the compressional wave velocity 
and attenuation. In addition, the receiver 2 may be translated 
along a line parallel to the sample surfaces, in order to 
capture the incident shear waves where their amplitudes are 
largest (see FIG. 3). The alignment angles of both trans- 
mitter and receiver are preferably chosen identical. The 
position of the receiver, denoted by x s is chosen according 
to a formula: 

X l S T tan 6+S R tan &+h tan 6, 

[0057] where 1 means either the shear (S) or compressional 
(P) waves, &J- and S R are respectively the transmitter and 
receiver stand-offs from the sample nearest surface (mea- 
sured when 6-0), and 8 their alignment angles. This formula 
is such that the axis of maximum amplitude of the shear 
beam transmitted through the sample impinges on the 
receiver aperture close to its centre (as per FIG. 3) because 
it is the place where the efficiency of said receiver is optimal. 
In this formula 8 S is the axis of the transmitted shear beam: 




[0058] Here V s and V f refer respectively to the velocity of 
the shear wave in the cement and sound speed in the medium 
in the absence of the sample. So for each incidence angle, 
the receiver is positioned at x s to capture the shear wave. 
Since the wave velocities in the cement are not known, one 
can use a guess value for the shear wave velocity and 
measure signals at various angles and receiver positions. The 
model-based method does not rely on accurate estimation of 
x s since it can account for arbitrary alignment and position 
of the receiver as described below. 

[0059] FIG. 4 represents an other embodiment of a mea- 
surement apparatus. In this embodiment, neither the trans- 
mitter and receiver nor the cement sample are rotated. 
Actually, the transmitter is constituted by an array of trans- 
ducers that each emits acoustic waves that are spaced of a 
short period of time from one transducer to an other. 
Consequently, the wave front that is emitted towards the 
sample is globally inclined of an angle 8. This apparatus also 
has electronic means that allow to change the period of time 
spacing each wave, in order to modify the incidence angle 
8 of the front wave. In an other embodiment, it is the 
receiver that is constituted by an array of transducers that 
each detects acoustics waves that are time spaced from one 
transducer to another. 

[0060] The embodiments described above refer to part of 
the ultrasonic subsystems. Accompanying these subsystems 
are ultrasonic pulse r (frequency about 500 kHz), signal 
conditioning circuitry including low and high pass electronic 
filters, amplifiers, as well as mechanical subsystems capable 
of rotation and/or linear translation of the sample and/or the 
transducers when required. Additionally, special transducer 
housings or windows may be necessary to shield them from 
the high temperature and pressure measurement environ- 
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meat. If such devices are used, the theoretical model 
described here above may be extended to account for 
transmission of the acoustic waves through them 

[0061] After being detected by the receiver, the signals are 
compared to calculated signals provided by the calculating 
method. At the end of the iteration process, the method will 
extract the velocities and attenuations of the compressional 
and shear waves from ultrasonic signals acquired in a 
cement sample and finally provide the mechanical charac- 
teristics of said sample. The velocities are denoted by Vp 
and Vs and the attenuations by ap and as. The calculating 
method is based on the use of a theoretical model that 
simulates the measurement technique implemented in the 
apparatus in order to provide signals corresponding to the 
measured signals and to compare these two sets of signals. 

[0062] The inputs to the model consist of several param- 
eters that describe the components of this measurement 
apparatus such as transducers, cement sample, medium fluid 
that separates transducers from the sample, as well as the 
system electronics. Each of these parameters is estimated. 
The velocities and the attenuation properties of the sample 
of tested material also constitute inputs parameters. These 
are estimated and initialised in the calculating method. This 
theoretical model is implemented in a computer code using 
known programming languages. Additionally, an interface is 
developed, said interface permitting to drive the calculating 
code automatically and with minimal operator intervention. 

[0063] For each set of measured signals, the computer 
code generates corresponding calculated signals. Certain 
attributes of the two sets of signals, such as transit times and 
amplitudes of particular events, are compared and differ- 
ences noted. On the basis of these differences, the code 
specifies updates in the input parameters (one or few or all 
at a time) as detailed below. The new set of inputs is used to 
generate a second set of calculated signals and comparisons 
in the attributes are performed again, and so on. The iteration 
process stops when the differences in all the attributes are 
within a tolerance factor prescribed by the operator or when 
the number of iterations exceeds a prescribed maximum. 

[0064] Optimal utilisation of the theoretical model and 
numerical code to calculate the cement parameters requires 
two calibration steps: 1) a characterisation of the radiation 
and reception properties of the transducers and 2) estimation 
of the transfer function of the electronics system which is not 
modelled, but is accounted for in the method according to 
the invention. 

[0065] The transducer characterisation may be performed 
according to well known techniques in a laboratory setting 
using a water tank. In this measurement, the receiver scans 
the acoustic field radiated by the transmitter. The scan is 
performed on a plane preferably perpendicular to the trans- 
ducers sensitivity lines. Then, the time-domain voltage 
signals are Fourier transformed to the frequency domain so 
as to obtain the transducers characteristics at each frequency 
component of the meaningful spectral range (which typi- 
cally spans 50 kHz to 1000 kHz). The frequency-domain 
voltage serves as reference data in a mathematical model 
characterising the measurement of the transducer. By fitting 
the calculated data to the measured data at each frequency, 
the optimal width and distribution of the transducer aper- 
tures are obtained. These values constitute inputs to the 
calculating model for the cement properties. 



[0066] Since environmental parameters such as tempera- 
ture and pressure influence the responses of the transducers 
and their electronics, their transfer function is captured from 
signals gathered preferably with the cement measurement 
apparatus. For this purpose, the measurement of the propa- 
gation of the signal in the fluid medium without the cement 
sample is used. The normalisation step that permits extrac- 
tion of the transfer function requires dividing the temporal 
spectrum of the measured signal by that of the calculated 
signal. An additional result of this calibration step is the 
estimation of the sound speed in the medium fluid. If the 
cement measurement is performed at various temperatures 
and pressures, then this calibration step (capturing of the 
transfer function) is required for each different set of tem- 
perature and pressure. Such a task may be accomplished 
without removing the sample from the apparatus but rather 
pushing the sample away from the line of sight of the 
transducers such that there is a clear fluid path for the signal, 

[0067] The cement characterisation model still needs two 
additional inputs (besides the four parameters that have to be 
estimated (Vp, Vs, ap, and as)). These are the cement 
density and its thickness (measured along the direction of 
signal propagation). Both of these parameters are measured 
independently using well known techniques such as buoy- 
ancy method for the density. If it is estimated that variations 
in temperature within the measurement apparatus leads to 
variations in the cement thickness that are larger than the 
estimated errors on its measurement, then this variations 
should be measured and taken into account in the calculation 
of the cement parameters. 

[0068] After this calibration step is performed, the math- 
ematical model according to the invention can be used. Once 
the computing system has provided — at the last iteration — 
the values for the velocities and attenuations for both shear 
and compressional waves signals from the receiver, these 
values have to be converted into known mechanical prop- 
erties to determine the integrity of the cement. The math- 
ematical relations between velocities and mechanical prop- 
erties are as follows: 

[0069] Let \ and pt be respectively the well-known Lame 
constants (// is also known as the shear modulus) of an 
elastic solid and E, v, and k, be Young's modulus, Poisson's 
ratio, and bulk modulus, respectively. These quantities are 
related to Vp (velocity of compressional waves), Vs (veloc- 
ity of shear waves) and density p of the tested material 
through the following: 

A=p-(V p 2 -2*V, 2 ) n=p-vh PI 
/i-(3A + 2>i) A t , 2 



[0070] In the data reported by the model, E, /* and k are 
expressed in units of pounds per square inch (psi). This is 
related to the SI unit Pa [kg/(ms2)] via the following 
relation: 
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ElPa) 
E(PSi) = 6.894757. 103 



[0071] Now, the mathematical model used to implement 
inputs in the method according to the invention will be 
described in more detail: 

[0072] In a preferred embodiment, the experiments take 
places inside a fluid medium (water for example) and the 
transducers are ultrasonically-pulsed acoustic piston ones 
propagating acoustic waves in transmission mode in a finite 
planar sample of cement (of thickness preferably between 10 
mm and 20 mm). 

[0073] Examples of a possible configuration and of its 
geometrical version where the cement layer is laterally 
infinite are shown in FIGS. 5(a) and (b). A Cartesian 
coordinate system (x, y, z) is used with the y axis pointing 
out of the plane of the paper. For this description of the 
model, the cement layer is assumed to be isotropic and 
homogeneous with propagating waves in the linear regime. 
However, these restrictions do not constitute a limitation of 
the calculating model as this latter can accommodate more 
advanced mathematical models for wave propagation in the 
cement layer. The transducer apertures have finite sizes in 
the (x, y, z). The transmitting (T) and receiving (R) trans- 
ducers may be modelled via the known Complex-Trans- 
ducer-Point (CTP) technique described for example in "A 
complex-transducer-point model for emitting and receiving 
ultrasonic transducers" from S. Zeroug, F. E. Stanke and R. 
Burridge in Wave motion Vol. 24, 21-40 1996. The CTP 
technique simulates piston transducers with Gaussian par- 
ticle velocity distributions. The same technique can also be 
generalised to accommodate transducers with arbitrary dis- 
tribution by weighted sums of these CTPs (ie, transducers 
modeled with the CTP technique). However, the formula 
given below [5] is general to account for arbitrary transducer 
aperture distributions and not restricted to CTPs. 

[0074] The forward mathematical model calculates the 
receiver voltage due to the pressure wave transmitted 
through the cement. For this purpose, spectral plane-wave 
expansion and synthesis of all pertinent wave fields are used 
to solve for the radiation, propagation in the fluid and 
transmission through the cement, and reception processes. 
The spectral integral representing the receiver voltage is 
computed via a real-axis integration scheme. In the follow- 
ing, the main formula of the model are presented. 

[0075] The time-domain voltage e^t) generated by the 
transducer (R) in reception mode can be derived from its 
frequency-domain counterpart E R (oo) as follows: 



e R (0= —J ERiuye-^du 



[0076] where (o is the angular frequency (w=2nf, where f 
is the frequency). 

[0077] In the calculating model, equation [4] is carried out 
via a Fast Fourier Transform (FFT) algorithm E(co) is given 
by: 



[0078] where 



-r 



[0079] and 



[0080] are the spectral plane-wave amplitudes of the nor- 
mal components of the particle velocity wavefield across the 
transducer aperture, T(k) is the spectral plane wave trans- 
mission coefficient that accounts for the acoustic wave 
interaction with the cement layer, and <!> is the "propagator" 
factor that translates the signal from the transmitter aperture 
to the sample entry surface and from the sample exit surface 
to the receiver aperture. In equation [5], yto) is a frequency- 
dependent quantity introduced to account for the transfer 
function of the transmitter and receiver electronics; it is 
captured through a calibration procedure as described above. 
The propagator O is given by the relation: 



2 

c expj -*r) + - jy) + -h - Z T )]} 



[0081] Various methods are available to express 



[0082] for a rectangular or apodised velocity distribution 
of a piston transducer or a phased array with radiation or 
reception characteristics close to that of a piston transducer. 
In a preferred embodiment, radiating (T) and receiving (R) 
CTPs are used. In this case, E R (co) is conveniently expressed 
in the following equation: 



with: 
k 3 =fcj + ij; 
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-continued 

ID 



[0083] where p f is the fluid density, Vf is the sound speed 
in the fluid, h is the cement layer thickness. The quantities 
(Xt r , Y t r , 2r r) are the CTP coordinates of T and R. [A 
detailed analysis can be found in the reference by Zeroug et 
al which is cited above.] For fixed transducers located and 
aligned in the y*0 plane as shown in FIG. 5 (b) with the 
sample rotated at angle 9, these complex coordinates can be 
written as: 



y T = yr = 0 y R = y* = 0; 



2r = - 
with 



dzv%9-h 



dcos8 - h 



- ibgcos8 



[9] 



[10] 



[0084] Here, W T R is the waist 1/e radius of the rotationally 
symmetric 3D Gaussian beams associated with T and R. 
From equation [11], bx and b R are identified as the Fresnel 
lengths of T and R. Since the real piston transducers have 
identical geometrical dimensions, we will assume that 
bT=b R . 

[0085] The expression for T(kx, ky) for an isotropic and 
homogeneous solid medium is given by: 



iY(A + S) 
(A+iYXS-iY) 



[12] 



[0086] where: 

A-Aty)-^-!*?) 2 tan(K p *#)+4jfc 2 K,K p tan(K,/i/2) 

[0087] 



[13] 



ric modes in the cement layer. The attenuations may further 
be expressed as linearly dependent on the frequency: a p> 
s-(u>/2tt)-o^ s ; with a slope given by and a, (expressed in 
Nepers/m/Hz). This function is preferred because it easily 
fits the behaviour of conventional cements. However, other 
frequency dependencies, perhaps needed for certain special 
cements, can also be accommodated in the calculating 
method according to the invention. 

[0089] The mathematical equations above are imple- 
mented numerically to provide the calculating model, which 
is used to generate signals corresponding to the measured 
signals. 

[0090] In practical, the acoustic method according to the 
invention proceeds as the following: 

[0091] Attributes are estimated from both the model - 
calculated and measured signals, said signals being voltages 
values as a function of time. In the preferred embodiment 
according to the invention, these attributes are: 

[0092] transit times of the P (compressional) and S 
(shear) events which may be estimated from the first 
zero crossing. 

[0093] amplitudes of the P and S events which may 
be defined as the peak-to-peak amplitude, the peak 
amplitude of the envelope of the analytical signal 
associated with the event, or the area underneath this 
envelope. 

[0094] The attributes for the calculated and measured 
signals are compared and the differences used to update for 
the model calculations the values of the velocities and 
attenuations for the acoustic waves trough the sample of 
material. This iteration process stops when the differences 
reach a prescribed tolerance factor or a maximum iteration 
number. Those two possible stopping keys are previously 
chosen and implemented in the calculating model. In a 
preferred example, simple linear minimise rs that have been 
found effective and robust for the various tests conducted are 
used. Other minimisation schemes such as based on least- 
square may be used. 

[0095] For the velocities, the calculating iteration is the 
following: let, t, ref be the transit time of the compressional 
waves event (P), in which case 1«P or of the shear waves 
event (S), in which case 1«S, in the measured signals. Let, t„ 
be the corresponding transit time of the calculated signals 
with a certain value for the wave velocity, V,. At the i th 
iteration (i=N)» Vj is updated according to: 



A = A(k)= [k] -2A 2 ) 2 tan(* P A/2) + 4k 2 K s K p tan{*,h/2) 
S = S{k) = {k] - 2k 2 ) 2 col(K p h/2) + ^K.K^o^hjl) 

Y=Y{k) = pk* s K p {K f 

and 



k PA = uj/ Vpj + i&pj 



[13] 



[14] 



[0088] Here Vp, Vs, dp,^ are the compressional and shear 
wave velocities and attenuations (expressed in Nepers/m), 
respectively, inside the cement, p is the fluid-to-cement 
density ratio pf/p c , and S and A are known as the Lamb 
characteristic functions for the symmetric and anti-symmet- 



v; +1 =v; + ,Sv ( ; 

where 
<5v,=- 



[15] 



[0096] where 1 means either the shear (S) or compressional 
(P) waves, d is the inter-transducers distance, h is the cement 
sample thickness, and V f is the sound speed in the fluid, all 
assumed known from independent measurement or calibra- 
tion steps. The iteration process stops when the relative 
difference between two succ 
prescribed tolerance factor: 
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[0097] The iteration process also stops if the number of 
iterations exceeds a prescribed maximum. 

[0098] For the attenuations, the calculating system is the 
following: let Aj rcf be the amplitude of the compressional or 
shear waves events in the measured signals (as per the 
definition given above. Let A! be the corresponding ampli- 
tude of the calculated signal with a certain value for the 
attenuation a,. At the i th iteration (i^N) <x lf is updated 
according to: 

oTooJ + d^; in] 

where 




operation of a transmitting element to a receiving element or 
that of a transmitting array to a receiving array under the 
assumption that the beam of the transducer array has similar 
propagation properties to the beam radiated by a single 
transducer of width equal to several wavelength (the refer- 
ence wavelength is measured in the fluid medium at the 
centre frequency of the transducer signal). 

[0103] Different tests have been prosecuted at ambient 
temperature and atmospheric pressure on solid samples. 
Several cements have been tested: 

[0104] 1 . Classical cement (API class G): cement and 
water; 

[0105] 2. lightweight cements at various densities: 
formulations including micro-cement and hollow 
microspheres of alumino silicate and water; 

[0106] 3. Foam cements at various density: classical 
cement+foaming agent and air 

[0107] 4. Composite cements (flexible cement): 
cement and flexible particles (XE 900 and XE 911) 





Cement 






Attp 


Atts 








Density 


Vp 


Vs 


(dB/cm/ (dB/cm/ 


E 




Description 


(kg/m 3 ) 


(m/s) 


(m/s) 


MHz) 


MHz) 


QqsQ 


V 


Classical cement (API class G) 


2025 


3736 


2040 


1.24 


5.36 


3309 


0.29 


Flexible cement with XE901 


1646 


2923 


3558 


2.99 


6.73 


1509 


0.3 


Flexible cement with XE900 


1633 


2538 


3293 


4.68 


33.33 


1046 


0.32 


Foam cement (30% air) 


1388 


2248 


3287 


32.93 


32.68 


838 


0.26 


Foam cement (20% air) 


1420 


2393 


3379 


8.3 


11.98 


743 


0.3 


Lightweight cement 1 


1819 


3358 


3797 


0.37 


7.69 


2347 


0.26 


Lightweight cement 2 


1595 


3064 


3694 


3.4 


7.77 


3700 


0.28 


Lightweight cement 3 


1470 


2848 


3496 


0.89 


4.45 


3249 


0.33 


Lightweight cement 4 


1116 


2786 


3429 


0.68 


6.52 


873 


0.32 



[0099] where 1 means either the shear (S) or compressional 
(P) waves. 

[0100] The iteration over the attenuation stops when the 
relative difference between two successive iterates of a 3 is 
within a prescribed tolerance factor: 




[18] 



[0101] The iteration process also stops if the number of 
iterations exceeds a prescribed maximum. 

[0102] The mathematical model has been described for 
generating calculated signals corresponding to a single 
transmitter and a single receiver interacting with a cement 
sample immersed in a fluid medium under general condi- 
tions for the transducer and cement sample alignment angles 
and positions. In case of a measurement with an array of 
transducers (as previously described with FIG. 4), this 
model is adequate to calculate signals corresponding to the 



[0108] These show that the method according to the 
invention simply provide accurate mechanical property val- 
ues of various material. Said mechanical properties can then 
be used as inputs in mathematical models to interpret data 
from known acoustic tools. 

[0109] Different examples have also been handled with a 
set of measured signals acquired from a sample made of a 
conventional class H (API system) cement. FIGS. 7-12 show 
comparison, along different iteration steps, between calcu- 
lated signals according to the method of the invention and 
measured signals. 

[0110] FIG. 6 shows the entire set of signals acquired with 
the measurement apparatus. The lower signal in the figure, 
denoted by "Fluid path", corresponds to the absence of the 
cement sample. This signal is used for calibration purposes 
of the method (as previously explained). The rest of the 
signals are acquired in presence of the cement sample when 
rotated at various angles ranging from 0° to 35°. The signal 
at 0° may be used to extract the compressional wave velocity 
and attenuation. Signals at 20° and above feature shear wave 
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events. One or more than one of these signals may be used 
to calculate the shear wave velocity and attenuation. 

[0111] FIGS. 7, 8 and 9 show comparisons between the 
measured and model-calculated signals of compressional 
wave velocity and attenuation for a sample rotation of 0° 
generated at the first, second, and final iterations of the 
calculating algorithm The tolerance criterion for both veloc- 
ity and attenuation is selected at 1%. The optimal values for 
the final compressional wave properties of the cement are 
written in the title of the plot in FIG. 9. 

[0112] FIGS. 10, 11, 12 show comparisons between the 
measured and model-calculated signals of shear wave veloc- 
ity and attenuation for a sample rotation of 25° generated at 
the first and final iterations of the calculating algorithm for 
the shear wave parameters. It can be noted that the discrep- 
ancy in the early compressional wave events is due to a 
limitation of the transducer characterisation model, which at 
this sample rotation affects only the compressional wave 
events and not the shear wave events. Also, the absence of 
the compressional wave event in the calculated signal at the 
first iteration is done purposely for shear event processing 
convenience. The tolerance criterion for both shear velocity 
and attenuation is selected at 1%. The optimal values for the 
final acoustic properties of the cement are written in the title 
of the plot in FIG. 12, the mechanical properties such as the 
Young's modulus, shear modules and Poisson's ratio are 
also calculated according to these final iterated acoustic 
characteristics. 

[0113] Once the four mechanical parameters are calcu- 
lated by the method for on the basis of two signals (a single 
rotational angle), a quality check can be performed using 
comparisons for signals at other sample rotation angles. 
FIG. 12 displays the results of such a quality control 
process. 

[0114] The method according to the invention thus pro- 
vides, from acoustic measurements through the cement 
sample, precious indications on the integrity of said cement, 
its thickening time in various temperature and pressure 
conditions. This method is consequently adapted to provide 
assessments on whether a particular cement is suitable for a 
particular cementing operation or how a particular cement 
will behave under particular well conditions. 



We claim: 

1. A method for estimating the time varying mechanical 
properties of a material comprising: 

propagating acoustic waves through a sample (3) of said 
material; 

measuring signals corresponding to said acoustic waves 
after they have propagated in said sample (3); and 

calculating the time varying mechanical properties of said 
material; the method characterised by: 

comparing attributes from said measured signals with 
attributes from corresponding signals provided by a 
calculating model, said calculating model account- 
ing for at least initially estimated acoustic properties 
of said material; 



using the differences between said measured and cal- 
culated signals to update said initially estimated 
acoustic properties of said material in the calculating 
model; 

processing (N) iterations of said comparisons and 
updatings until a given calculating stop; and 

calculating the time varying mechanical properties of 
said material from said final iterated acoustic prop- 
erties. 

2. A method according to claim 1, wherein the acoustic 
waves are both compressional and shear waves. 

3. A method according to claim 1, wherein the given final 
stop is when the differences between measured and calcu- 
lated signals are within a given tolerance factor. 

4. A method according to claim 1, wherein the given final 
stop is when the iteration number (N) reaches a prescribed 
maximum. 

5. A method according to claim 1, wherein the attributes 
of the measured and calculated signals are the transit times 
of the acoustic waves and amplitudes of said signals. 

6. A method according to claim 5, wherein the acoustic 
properties of the material comprise the velocity of the 
acoustic waves through the sample (3) of said material. 

7. A method according to claim 1, wherein the propaga- 
tion and measuring means comprise ultrasonic compres- 
sional transducers, among which a transmitter (1) and a 
receiver (2). 

8. A method according to claim 7, further comprising the 
steps of immersing the transmitter (1) and the receiver (2) in 
a medium (4) and placing the sample (3) of material between 
said transmitter and said receiver. 

9. A method according to claim 8, wherein at the I th 
iteration (where i^N), the value of the velocity for the 
acoustic waves through the sample (3) of material is updated 
according to: 



Vr = v; + 6v ( \ where 6 V( = V/ 

where 1 means either the shear (S) or compressional (P) 
waves, \{ ct is the transit time of the waves in the measured 
signals, t, is the corresponding transit time of the calculated 
signal with a certain value V, for the wave velocity, d is the 
distance between the propagating and the measuring means, 
h is the sample thickness, and V f is the sound speed in the 
medium. 

10. A method according to claim 9, wherein the iteration 
stops when the relative difference between two successive 
iterated values of the velocity for the acoustic waves through 
the sample (3) of material is within a prescribed tolerance 
factor e V[ such as: 

\v,\ <ev >- 

11. A method according to claim 1, wherein the time 
varying mechanical properties of the material are the 
Young's modulus, the shear modulus and the Poisson's ratio. 
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12. A method according to claim 1, wherein the acoustic 
properties of the material comprise attenuation of the acous- 
tic waves through the sample of material. 

13. A method according to claim 12, wherein at the I th 
iteration (i§N) 5 the value of the attenuation for the acoustic 
wave is updated according to: 

.. . 

of = aj + c5o.; where 4», = *i 

A? 

where 1 means either the shear (S) or compressional (P) 
waves, Ai ref is the amplitude of the acoustic waves in the 
measured signals and A 1 is the corresponding amplitude of 
the calculated signal with a certain value for the attenuation 
<V 

14. A method according to claim 13, wherein the iteration 
stops when the relative difference between two successive 
iterated values of the attenuation for the acoustic wave is 
within a prescribed tolerance factor e Vl a such as: 



15. A method according to claim 1, wherein the calculat- 
ing model accounts for characteristics of the propagating 
and measuring means. 

16. A method according to claim 15, wherein the charac- 
teristics of the propagating and measuring means comprise 
the radiation and the reception properties of said propagating 
and measuring means and the estimation of the transfer 
function of the electronic devices included in said means. 

17. A method according to claim 1, wherein the calculat- 
ing model accounts for environmental parameters such as 
temperature and pressure. 

18. A method according to claim 18, wherein the material 
is a settable material. 

19. An apparatus for propagating and measuring acoustic 
waves through a sample (3) of material, the mechanical 
properties of which vary with time, said apparatus compris- 
ing: 

at least one transmitter (1) propagating said acoustic 
waves towards at least one receiver (2), said receiver 
detecting said acoustic waves after they have passed 



through said sample; and means for calculating prop- 
erties of the sample based on the detected acoustic 
waves; characterised in that the apparatus comprises: 

means for varying the incidence angle (8) of the acoustic 
waves on the sample of the material; and 

means for 

comparing attributes from said detected signals with 
attributes from corresponding signals provided by a 
calculating model, said calculating model account- 
ing for at least initially estimated acoustic properties 
of said material; 

using the differences between said measured and cal- 
culated signals to update said initially estimated 
acoustic properties of said material in the calculating 
model; 

processing (N) iterations of said comparisons and 
updatings until a given calculating stop; and 

calculating the time varying mechanical properties of 
said material from said final iterated acoustic prop- 
erties. 

20. An apparatus according to claim 19, wherein the 
sample of material is rotational such that its surface makes 
different incidence angles (0) with respect to the transmitter 
and receiver sensitivity line. 

21. An apparatus according to claim 20, wherein the 
transmitter (1) and receiver (2) alignment is controlled via 
an external mechanism. 

22. An apparatus according to claim 21, wherein the 
receiver translates along a line parallel to the surface of the 
sample of material. 

23. An apparatus according to claim 20, further compris- 
ing multiple pairs of transmitters and receivers. 

24. An apparatus according to claim 23, wherein the 
transmitters and receivers are aligned from each part of the 
sample of material such that each transmitter-receiver pair 
sensitivity line makes a different incidence angle (8) with 
respect to said sample. 

25. An apparatus according to claim 23, wherein the 
number of pairs is such that several incidence angles (0) 
from zero to a maximum angle are covered, said maximum 
angle being preferably superior to 30°. 

26. An apparatus according to claim 19, wherein the 
transmitter or the receiver comprises an array of transducers 
that respectively each emits or detects acoustic waves that 
are time spaced from one transducer to another. 



* * * * * 
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